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Abstract. We consider a simple atomic two-body bound state system that is overall charge neutral and 
placed in a static electric and magnetic field, and calculate the magneto-electric response function as a 
function of frequency. This is done from first principles using a two-particle Hamiltonian for both an 
harmonic oscillator and Coulomb binding potential. In the high frequency limit, the response function falls 
off as l/oj^ whilst at low frequencies it tends to a constant value. 



PACS. 31.15.A- 
erties 



Ab initio calculations - 31.15.ap Polarizabilities and other atomic and molecular prop- 



1 Introduction 



Magneto-electric effects are by now a well established phe- 
nomena both theoretically and experimentally. The ten- 
dency has been to focus on relatively large molecular sys- 
tems where DFT calculations [TJ[2] apply and experimen- 
tal values have been measured [3]. In addition to being 
associated with important optical phenomena [2113], their 
existence has played an important role in the Casimir 
physics [S1[SJ[7], in particular if it is possible in certain 
circumstances to find a contribution from the quantum 
vacuum to a bodies momentum. For these reasons we con- 
sider it an interesting question to ask what are the simplest 
models that describe possible and display magneto-electric 
effects. 

In this article we calculate the magneto-electric re- 
sponse function for the two simplest bound state systems - 
the harmonic oscillator and the hydrogen atom. The prin- 
cipal difference between these two atomic systems is that 
the harmonic oscillator is strongly bound whereas the hy- 
drogen atom with its Coulomb potential is weakly bound. 
This manifests itself in the accessibility of different energy 
eigenstates in the perturbation theory 

This article is organised as follows. In Section [5] we de- 
fine the atomic system and formulate the response for an 
arbitrary binding potential. As a tractable example that 
has a closed form, the harmonic oscillator is chosen and 
its magneto-electric response function calculated explic- 
itly. In Section [3] the response function for hydrogen is 
presented. It is then compared to the appropriate DFT 
result and its relation to experimental values. Finally, in 



Section|3]we summarise our results and provide some com- 
ment on their validity and applicability. 

2 General formulation 

We will now derive the ME response function for a two- 
body charge neutral composite system in static external 
fields (E", B°) using the Hamiltonian formalism (the bold- 
face used here is to indicate they are external fields) . One 
can also use the path integral approach and the coupled 
classical Lorentz force equations to obtain information 
about the response function. The path integral gives a cor- 
rect response function at high frequencies but it is not reli- 
able at low frequencies as well as suffering from the wrong 
analytic structure in the complex frequency plane jSj. A 
consideration of the coupled Lorentz force equations pro- 
duces similar difficulties and so we use the Hamiltonian 
method exclusively in this paper. 

The system we consider is illustrated in Figure [1] Two 
equal but opposite electrical charges with coordinates {q\, ^2) 
and masses (toi, 7712) interacting with a classical (c-number) 
gauge field (</>, A') that are the combined contribution of 
the static external fields and the fluctuating source that is 
used to probe the system. A binding potential V{qi — (72) 
holds the charges together giving a bound state that is 
overall charge neutral. The Hamiltonian that describes 
this system is given by 



H 



1 



2toi 
1 



{pi ~ eA{q,)f + e(t>iqi) 



+ o + e^fe)) - 60(92) + V{qi - q2). (1) 

2to2 
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[m^,-e,q^ 




Fig. 1. The bound state two-body atomic system consisting 
of two equal but opposite electrical charges, with two different 
masses mi and m2 and coordinates q\ and ql, (with respect to 
the coordinate system E). To deduce the response function for 
this system we will change coordinates to the centre of mass 
coordinate X' and the separation coordinate a:'. 



To calculate the response function of this body due to a 
high frequency electromagnetic field, it is necessary to pass 
to a new set of coordinates that consists of both the centre 
of mass coordinate X"^ and the separation vector x*. These 
are shown in Figure [TJ The new variables are then defined 
by X := {miqi+m2q2)/M, x := {qi -172), Af := rrii -f m2, 
m :— 11111712/ M , and niA = — mi. Correspondingly, we 
change from the two particles momenta (pi,p2) to {P,p) 
where P is the conjugate momenta of the centre of mass 
coordinate X and p is likewise the conjugate momenta to 
the separation vector x. It is also necessary to implement 
this change of coordinates on the gauge field; it can then 
be expanded about the centre of mass coordinate as 



A{X + {m2/M)x) 



A{X)+[^)x^V.A{X) 



(m2\' 



2) \ M J 
mi 



'x^V,V,A{X), (2) 



AiX - {mi/M)x) = A{X) - l^^j x^V^iX) 

I) (^)'^*^'V,V,A(X), (3) 



and similarly for the scalar potential cf). We will work to 
first order in the spatial derivatives of the gauge potential 
and thereby neglect the last terms in the above expansion. 
This corresponds to electric and magnetic fields that can 
vary in time, but that are spatially constant (so the ap- 
proach is restricted to wavevectors that are less than the 
inverse of the size of the atomic system) . 

Using the Lagrangian as an intermediate step in per- 
forming the change of coordinates, we recognize here that 
E{t, X) = -V(f>{t, X) - dtA{t, X), whilst the derivative of 
the vector potential once projected with the Levi-Civita 
tensor will give the magnetic field. The Hamiltonian in the 



new coordinates after making this expansion reads 
2m 2m 



+ ^ ["^y x'x^V^Ak{t,X)V,Akit,X) 
+0{V'^A). 



(4) 



The ME activity results from a source magnetic field in- 
ducing electrical polarisation. We can now define the magneto- 
electric response of the bound state system by promoting 
all of the canonical degrees of freedom to operators. It is 
defined by 



{e6x,{t)) 



dt'xtr{t~t')SB,{t',X), 



(5) 



where the Ihs is the standard expectation value of the 
fluctuating electric dipole moment induced on the rhs by 
an externally applied fluctuating magnetic field (i.e. a test 
source) that is in general time and space dependent. In 
terms of correlation functions, it is given by the retarded 
two-point function 



{t-t') = -ie{t - t'){Q\[e5x,{t),5dj{t')]\Q), (6) 



where the operator Ojit') couples to the fluctuating mag- 
netic flcld 6Bj{t' , X), and is to be found from the micro- 
scopic theory given by the Hamiltonian Equation (0]) . The 
ground state \n) that we will use will be specifled later 
in this section. Taking the Fourier transform of this one 
flnds in frequency space 



H - En-huj 



-(f2\exi 



1 



H - Eo + huj 



d,\ny 



(7) 



where Eq is the ground state energy of the system. To ob- 
tain the detailed form of the response function, we need 
to specify the form of the operator Oj . One can see from 
Equation (|4|) that the operator has two contributions. One 
is linear and the other is quadratic in the gauge poten- 
tial. To calculate the operator we make the specific gauge 
choice for the gauge potential 



-SB{t))AX. 



(8) 



With this choice the Hamiltonian is then a function of only 
gauge invariant quantities. As a final step to fully specify 
the Hamiltonian we include the static electric field by sub- 
stituting Ei{t, X) = E°. The final form of the Hamiltonian 
that will be used for subsequent calculations is given by 



H 



1 

2M' 



1 
2m 



P 



V{x) - ex' ■ E° 
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-e (^) (iAP).(B" + 5i3(t)) 



2Mm 



L-{B° + SB{t)) 



^im^jn]: (9) 



where L — x A p is the orbital angular momentum about 
the centre of mass origin. One can now just read off the 
operator that couples to the fluctuating magnetic field 



2mM 



L, 



4to V~M~/ 



X°'x''[SabSij 



(5„<5fcj]B°. (10) 



There is also a contribution from the term linear in the 
centre of mass momenta. However, this is zero when eval- 
uated between the ground states (zero momentum plane- 
wave eigenfunctions) . This amounts to a choice of refer- 
ence frame (the centre of mass frame) and can be used to 
define the components of the static external electromag- 
netic fields as well. Indeed, we have not specified so far the 
nature of the ground state | f2) . For computational conve- 
nience we work with a perturbed ground state due to the 
presence of static external electric field 



\f2) = |o) -ce: 



\ ' 1 

n—1 



\n){n\x'\0). (11) 



This choice corresponds to the physical situation where we 
put the two particle system in the static external electric 
field first, let the system settle down and then apply the 
static external magnetic field. From Equation pUj) we see 
there are two separate contributions to the response func- 
tion. The first is due to a coupling with the angular mo- 
mentum operator which we write as xfj^ {^i L) \ the second 
is due to a quadrupole moment coupling which we write as 
xfj^{i^,x^)- Then the full response function is just given 

by their sum x^^(w) = x^^(w,^) + xfj^(w,i^). First 
consider evaluating the contribution due to the angular 
momentum operator 



e Am 
mM 



{n\x. 



1 



mM 



■ {n\a 



H -Eo-hu 
1 

' H - Eo + hu} 



Lj\n) 

L,\QY.{12) 



X {Q\x- 



1 



1 



Hq — Eq + hoj Hq — Eq + huj 

Inserting two complete sets of states gives 



4|i7)*(13) 



3^2 



\Am^M^ 



{[n\xi\m){m\Lk\n){n\Lj\n) 

1 1 
^ 

E,n — Eq — flUJ En — Eq — flLO 

+ {{f2\xi\m){m\Lk\n) {n\Lj\Q))* 

:)■ 



1 



1 



Em — Eo + huj En — Eq + huj 

For the quadrupole moment contribution we have 

x{x^B° - B" • xxj)\n) 
/ e^m^ \ , / 1 



(14) 



, , , (f2\x. 

x{x^B"j -B'^ ■xxj)\n)*. (15) 
Inserting a single complete set of states gives 

3^2 



En — huj 



x(f2\xi\n)(n\{x'^B° - B° • xxj)\f2) 



)^{En-Eo 



V4mM2 J ^ \ Er,.- En + hu , 
x{{n\x^\n){n\{x^B°j - B° • xxj)\n))* . (16) 

When the perturbed ground state given by Equation ([TT|) 
is substituted into Equations (IT4|) and ((T6)) we arrive at 
a final expression for the magneto-electric response func- 
tion. To linear order in the static external fields Equa- 
tion (fT4ll becomes 



m.n s^O 



The next step is to expand the denominators in terms of 
the static magnetic field. From Equation ^ the perturba- 
tion of the Hamiltonian is given by SH — — (eTO/i/mM)LiB° 
therefore we find (keeping only the terms linear in the 
static magnetic field) 



X {il\xi 



1 



Hq — Eq — hui 



(ifeB°) 



nMJ 



mM j \ mM I 
1 



((0|x,|m)(m|Lfe|n)(n|Lj|s)(s|x'|0) 
1 1 1 

— Eq Em ~ Eq — hu! En — Eq — TlLO 

+mx^\m){m\Lk\n}{n\L,\s){s\x'\0))* 

1 1 1 

Eg — Eq Ejji — Eq -\~ hid Eji — Eq -\- Hlu 

The corresponding form that Equation fjl6p takes is 

4^2 



(17) 



Ho — Eq — Tllu 



/emA\ 



e mA 



mM I V mM J 



BjEf 



n S5^0 



Es — Eq En — Eq — hoj 
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x(({0\xi\n){n\{x 6kj — S:kXj)\s){s\x \0) coefficients. The full response function then takes the final 

+ {s\xt\n){n\{x'^Skj - XkXj)\0){0\x'\s)) 
mxi\n){n\{x'Sk, ~ Xkx,)\s){s\x'\0) xf.^'i^) = xf,"" i^, L) + xf,"" (co^x^) 



+ {s\x^\n){n\{x'Skj ^xkx,)\0){0\x'\s)y). (18) e*m^. , oj' 



2 I , ,2 



To go further, it is necessary to specify a binding po- 

tential so that the energy eigenvalues and eigenfunctions H — g 2 (-^i^i ^ I'^^'^i ^ I^-^" ' 

can be deduced. - a; 4 4 

(24) 



2.1 The harmonic oscillator binding potential 

As the simplest example one might consider, the harmonic 
oscillator with V{x) = mojQx'^ /2. This model is both pos- 
sible to solve analytically and relevant phenomenologically 
in displaying the properties associated with real matter. A 
notable feature here is that since the potential is strongly 
confining, the perturbed ground state due to the external 
static electric field has only the first excited state surviv- 
ing in the summation. Using the operator algebra of the 
oscillator and the defining relations 



i. = /^(al+a.) (19) 

Li = -ih/2eijkia'jak ~ a^al) (20) 
= -iheijka]ak, (21) 

the matrix elements can be evaluated explicitly in Equa- 
tion (ITTl) 



1 



1 



El — Eo — huj El — Eq — huj 
1 1 



El - Ea + fiuj El - Ea + hu 



4 9 

e 



(E'/B ' - (E" • B")%) 



(22) 



In an analogous use of the operator algebra the quadrupole 
contribution Equation (jl8p can be similarly evaluated 



1 



x(4E^B° - E^B" - (E° • B°)%). (23) 



Note here that both contributions have an anti-symmetric 
term upon writing the tensor structure in the external 
fields as a sum of symmetric and antisymmetric pieces, 
which will have a relevance to our later discussion. Both 
the angular momentum and quadrupolar contribution are 
of the same order as can be seen from their multiplicative 



3 The non-relativistic hydrogen atom 

We now consider the second simplest system, namely the 
hydrogen atom. This is a weakly bound system with a 
Coulomb potential given by V{r) — —e'^/{ATTeor) in spher- 
ical polar coordinates. Since the proton is very much more 
massive than the electron, we can make the replacement 
TO/i/M 1, leaving just the reduced mass in all ex- 
pressions (which is just the electron mass). Equation (1141) 
can be evaluated now using the hydrogenic eigenstates 
\n, L, m) and the energy spectrum En = Eijr? . One finds 



2 00 ria — l La 

B^EE E 

a— 1 ria— 1 La—Q m,a — — L 

1 1 

^ En^ — El — huj En2 — El — huj 

x((i7|xj|ni,ii,mi)(ni,ii,TOi|Lfe|n2,i2,TO2) 

x{n2,L2,m2\Lj\[2) 
1 1 



^Eni — El + Tvjj En2 — El + huj ^ 

((i7|£j|ni,Li,mi)(ni,ii,TOi|Lfe|n2,i2,TO2) 

x{n2,L2,m2\L,\f2))*. (25) 

For the quadrupole contribution Equation becomes 

Xtj''(.^,x') = +' 



Am? 



]b°. m^^\nl,Ll,ml) 



ni ,Li ,mi 



x{ni,Li,mi\S?5kj - XkXj\fi) 



1 



x{ni,Li,mi\x Skj - XkXj\f2)) 



Eni - El - huj 
{{^2\x^\nl,Ll,ml) 
1 



En^ — El + huj 



(26) 



Because the electron is weakly bound the perturbed ground 
state due to the external static electric field requires a full 
summation over the principal quantum number 



|r2) = |l,0,0)-eE?^ ^ 

n—2 m— — 1 

X (n, 1, m\xi\l, 0, 0), 



1 



En — i?l 



n, 1, m) 



(27) 
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in contrast to the strongly bound harmonic oscillator po- 
tential. In Equation (|25p the only non-zero elements are 
when Li ~ L2 — 1 (i.e. just standard selection rules or 
addition of angular momentum) and thus reduces to 



(-'^) = U^r"'^' E E E E 

^ ^ mi — — 1 7712 — — 1 n— 2 m— — 1 

1,0, 0|ii|n, 1, mi}(n, 1,to|x;|1, 0,0) (1, mi|ifc|l, TO2) 
X (1, m2\Lj\\, m)- 



1 



1 



(En -El- huY E„ - El 
+ ((1,0, 0|ii|n, 1, mi){n, 1,to|x;|1, 0,0) (1, mi\Lk\l,m2) 

x(l, m2|X, |l, m))*- -r . 

^ ' ' {En-Ei+hujfEn-Ei 

(28) 

Turning now to the quadrupole case. Equation ((26|) is 
slightly more complicated as the expectation values of the 
quadrupole operator requires performing some addition 
of angular momenta. Indeed, the quadrupole operator can 
be written as a linear superposition of spherical harmonics 
^2,771 and ^0,0- This implies non-trivial overlaps of matrix 
elements such that the summation over the L eigenvalue 
will not completely reduce to a single value as in the pre- 
vious case but rather one finds 



Am? 

E 



3 2 L 
,0 ^ 

ni—2 n—2 L—1 Li — 1 m- 



^l^^EEEE E 



1 



1 



E„ 



El En, -Ei-hu 



777 1 — — L 1 

x{{l,0,0\xi\ni, Li,mi){n, L,m\xi\l,0,0) 
x{ni, Li,mi\x'^6kj - XkXi\n,L,m) 
{l,0,0\xi\n,L,rn){n,L,m\xi\ni,Li, mi) 
x{ni,Li,mi\x'^Skj - ifeij|l,0,0)) 
1 1 



En -El En, -Ei+hu, 

x(((l,0, 0|ii|ni,Li,mi)(n, L, m|x;|l,0,0) 
x{ni,Li,mi\S?5k] - XkXi\n, L,m))* 
+ ((1,0, 0\xi\n, L, m){n, L, m\xi\ni, Li,mi) 
x{ni,Li,mi\x^6kj -xkx,\l,0,0))*). (29) 



Equations (|28p and ([29)1 are the main results of this pa- 
per. We see that like for the harmonic oscillator, there 
is a 1/w^ behaviour at large frequencies. To avoid the 
singularity when Tiuj = En — Ei, we supplement a phe- 
nomenological line width F so that any given excited en- 
ergy level can decay to a lower eigenstate by the replace- 
ment En — En + ir. We have evaluated Equations (^51) 
and analytically using Mathematica as an expansion 
in the principal quantum number. The real part of the sus- 
ceptibility is plotted in Figures [5] and O as a function of 
frequency both for the low and high frequency limits, and 
close to the resonance respectively. To do this we have used 
numerical data of typical field strengths of |B°| = lOT, 



Re(.Y[l,21)/(Vce„) 
7.x 10-'" 




■ ■ ' — - — ■ ■ ^ ' ' ' ^ ,, frequency (Hz) 

1.0x10" 1.5x10" 2.0x10" 

Fig. 2. The real part of the (dimensionless) off-diagonal 
susceptibility (the bi-anisotropic component X12 with V = 
(4/3)7rao the atomic volume in the ground state) of a hydrogen 
atom with parameters E" = (10^ l/m~\ 0, 0), B" = (0, 10r,0), 
bjQ = IQ^'' Hz and F = Kf Hz. This plot shows the zero and 
high frequency limits. 

Rc(^[1.2J)/(VcE„) 



10" 7.5x10" 7.5x10' 




frequency (Hz) 



Fig. 3. The real part of the (dimensionless) off-diagonal 
susceptibility (the bi-anisotropic component X12 with V = 
(4/3)7rao the atomic volume in the ground state) of a hydrogen 
atom with parameters E" = (10^1/m"\ 0, 0), B" = (0, 10r,0), 
luq — 10^^ Hz and F — IQ^Hz. This plot shows the resonant 
structure. 



|E^| = lO^T^m"^, a resonant frequency of wq = 10 
which corresponds to the energy difference between the 
first two levels in hydrogen, and a spontaneous decay rate 



of F 



10« 



In Figure 0] the imaginary part of the 



susceptibility is plotted close to the resonance. Of course, 
since the decay rate is very much smaller than the resonant 
frequencies there is a huge enhancement in the suscepti- 
bility close to resonance. 

It is worth pointing out here an issue of the validity of 
our calculation with respect to the photoelectric effect [9] . 
Given that we are driving an atomic system with an elec- 
tromagnetic wave at some frequency we may wonder if it 
is valid at high frequencies. In this regime the associated 
wave vector is also high and the photon can probe the 
shorter length scales and transfer more momentum to the 
electron. However, if one considers the differential cross 
section for photo-electric effect in hydrogen one knows at 
high frequency it behaves as . Therefore at high fre- 
quency the atom will remain intact and the previous anal- 
ysis of the magneto-electric response should remain valid. 
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Im(^[l,2])/(VCE„) 




Fig. 4. The imaginary part of the (dimensionless) off-diagonal 
susceptibility (the bi-anisotropic component xi2 with V — 
(4/3)7raQ the atomic volume in the ground state) of a hydrogen 
atom with parameters E° = 0), B° = (0, lOT, 0), 

uJo = W^^Hz and F = W^Hz. This plot shows the resonant 
structure. 

3.1 Size of the effect and comparisons 

It is instructive to give a numerical estimate of the size of 
the response function as compared to the standard electric 
susceptibility. To do this, consider the harmonic oscillator 
Equations and (1^ at zero frequency. Then up to nu- 
merical factors (i.e. simple dimensional analysis) we have 

~ (^-^) P (30) 

The last factor, /3, can be seen too be a dimensionless 
number, whilst the first is the standard static electric sus- 
ceptibility. Putting in the number leads to a factor of 
P ~ 10~^^. This will also be approximately the same for 
the hydrogen atom if the optical transition frequencies are 
chosen to coincide. 

We can compare this scale to experimental [3] and 
DFT values [T] found for the change in refractive An as 
compared to the absence of static external electromag- 
netic fields. For certain large complex molecules values 
of An = (iV/y)(x/(eoc)) ^ 10~^^ are found experimen- 
tally [3], where N/V is the number density of the sam- 
ple. For the helium atom a DFT calculation [T gives a 
refractive index difference of AnheUum ^ 10~^^ for a sam- 
ple with number density N/V ^ 10^^m~^ evaluated at a 
wavelength of A = 632. 8nm. The hydrogen estimate (for 
the same number density) from our calculations taking 
into account the number density scaling is Anhydrogen — 
{NlV){xl{toc)) - (7V/F)(eVeomcj2)/3 ^ iQ-is. It is also 
worth mentioning that |x^^(w)|, as for the harmonic oscil- 
lator, will scale as size of the system since it is proportional 
to the standard polarisability. For the helium atom, the 
static polarisability is a;jeiitnn(0) = 0.22 x 10~^°Coulomb 
meters/Volt, which when divided by eo gives a volume of 
16.6ag. For the hydrogen atom the corresponding volume 
is 4aQ from which we find a scaling factor of approximately 



four between the hydrogen and helium. The precise form 
however would have to be fitted empirically with the help 
of DFT calculations in order to match on to experimental 
values such as found in [3]. 



We have presented an exact quantum mechanical per- 
turbation theory calculation of the magneto-electric re- 
sponse function for atomic systems with the simplest bind- 
ing potentials, namely the harmonic oscillator and the 
Coulomb potential. We have deduced analytic forms for 
the magneto-electric response tensor as a function of fre- 
quency that can be calculated exactly. A common feature 
is that at high frequency they have l/w^ behaviour whilst 
at low frequency they tend to a constant value. It would 
be interesting to try and apply the same method exactly 
to the helium atom. 

There are interesting implications of this calculation. 
Concerning the Feigel effect .5 , where a net momentum 
density for a medium (Equation (21) in [5]) with an 
magneto-electric response function is developed, the re- 
sult found there is fourth power divergence in frequency, 
which is then simply cut-off. In fact it is only the anti- 
symmetric part of the magneto-electric response tensor 
that contributes to the momentum. This neglected how- 
ever the dependence on frequency and was treated as a 
constant. What we see now is that this divergence will be 
softened to a quadratic divergence though it will not be 
simply washed away altogether. This shows that the as- 
sumption of a cut-off at high frequencies for ME is not 
justified and that the divergence has to be resolved by 
other means, such as done recently in |10| . 
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